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Summary 

Since  the  overall  rate  of  fuel  evaporation  in  combustors  is 
strongly  influenced  by  the  gas: drop  relative  velocity,  detailed 
information  on  the  aerodynamics  of  the  flow  in  a given  combustor 
configuration  is  essential  for  its  prediction.  Such  aerodynamic 
information  can  be  obtained  either  from  measurements  or  from  mathematical 
simulation  of  the  events  taking  place  in  actual  systems.  Owing  to  the 
advances  made  in  the  field  of  mathematical  modelling  turbulent  flow 
phenomena,  and  the  development  of  powerful  techniques  for  solving  the 
partial  differential  equations  of  fluid  mechanics,  the  latter  approach 
is  on  the  waj^fco  becoming  moreeff icient  and  economical.  The  flow  equations 
can  be  solved  using  either  two  or  three  dimensional  models,  however,  almost 
all  practical  systems  may  be  treated  as  two-dimensional  (axisymmetric)  in 
the  important  region,  which  includes  the  region  of  soot  formation  near 
the  fuel  injector. 

The  calculation  procedure  is  here  applied  to  the  case  of  fuel 
injection  into  a recirculation  zone  behind  a baffle  in  a duct.  The 
elliptic  flow  equations  are  evaluated  using  the  k,E  model  of  turbulence. 

The  simultaneous  differential  equations  of  droplet  dynamics  are  evaluated 
in  each  cell  of  the  flowfield  for  each  size  group  of  the  spray  using  short 


time  steps.  To  adequately  represent  the  typical  fuel  spray  size  distri- 
butions, twenty  size  intervals  of  droplets  are  required. 

The  results  clearly  demonstrate  the  controlling  influence  of  the 
convective  contribution  to  spray  evaporation.  A particularly  interesting 
conclusion  of  the  analysis  is  that,  in  conditions  representative  of 
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combustion  systems,  the  droplets  have  just  deviated  from  their  initial 
conical  path  towards  the  flow  direction  when  they  vanish.  This  phenomena 
is  consistent  with  experimental  observations  and  to  a great  extent  justi- 
fies the  use  of  two-dimensional  modelling  for  the  spray  evaporation  region 
of  practical  combustors. 
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1.  Introduction 

At  the  present  time  about  one  half  of  the  World's  -fossil  fuel 
production  is  consumed  through  spray  combustion.  The  accurate 
prediction  of  the  combustion  behaviour  of  fuel  sprays  increases  the 
designers  ability  to  optimize  the  design  of  gas  turbine  and  industrial 
combustors,  diesel  engines,  etc.,  to  improve  efficiency  and  minimize 
pollutant  production. 

In  the  past  there  have  been  many  studies  on  the  vaporization  of 
single  droplets  (eg.  Ref.  1),  and  a smaller  number  of  studies  have  been 
concerned  with  the  combustion  of  clouds  of  droplets  (eg  Ref.  2) . In 
Refs  3 & 4 it  shown  that  in  most  practical  combustors  the  flame  does 
not  generally  stabilize  on  the  individual  droplets,  but  the  droplets 
first  evaporate,  then  the  fuel  mixes  with  the  air,  and  combustion 
actually  occurs  in  a turbulent  mixing  controlled  flame.  This  does 
not  mean  that  evaporation  does  not  play  a part  in  determining  the 
combustor  characteristics.  Indeed,  Fig  1 shows  that  the  measured 


stability  loops  for  a small  gas  turbine  combustor  vary  dramatically 
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with  the  fuel  nozzle  rating  i.e.  droplet  size  distribution  . This 
figure  clearly  demonstrates  that  the  finer  fuel  spray  (resulting  from 
a nozzle  with  a smaller  flow  rating)  gives  a greater  proportion  of  the 
fuel  evaporated  and  mixed  in  the  flame  stabilizing  primary  zone  of  the 
combustor. 

Recent  experimental  studies  have  shown  that  the  proportion  of 
pollutants  produced  by  gas  turbine  combustors  may  be  minimized  by  the 
use  of  pre-mixed,  pre-vapourized  fuel.  In  Ref  6 for  example,  it  is 
shown  that  pre-vapourizing  the  fuel  decreased  the  production  of  CO, 


NO^  and  smoke  number  by  30%,  24%  and  60%  respectively.  Pre-mixed 
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pre-vapourized  combustors  operating  with  lean  primary  zones  have 
achieved  NO  levels  down  to  0.3g  NO- /kg  fuel  at  the  supersonic  cruise 
condition^,  and  they  are  likely  to  be  suitable  for  future  synthetic 
fuels  having  a low  H/C  ratio.  A carburation  system  capable  of 
providing  a uniformly  pre-mixed,  pre-vapourized  fuel  is  described  in 
Ref  8.  There  is  thus  considerable  current  interest  in  predicting 
the  detailed  overall  behaviour  of  combustion  systems  with  various 
methods  of  fuel  preparation.  In  order  to  predict  quantitatively  the 
characteristics  (ie.  stability^ efficiency  and  pollutant  production) 
of  a combustor,  we  must  calculate  both  the  path  of  the  various 
drop  sizes  in  thfe  spray,  and  the  distribution  of  the  fuel  evaporated 
from  the  spray.  In  current  gas  turbine  combustors,  as  in  many  other 
applications  of  spray  combustion,  the  droplets  evaporate  along  their 
trajectories  in  che  presence  of  a forced  convective  motion  of  the 
gas,  leaving  behind  pockets  of  fuel  vapour  which  later  mix  with 
the  oxidizer  and  burn.  While  smaller  droplets  evaporate  almost 
instantaneously,  larger  ones  may  penetrate  the  flamefront  depending 
on  the  compromise  between  the  degree  of  atomization  of  the  liquid 
spray  and  its  issue  velocity.  Since  the  overall  rate  of  evaporation 
is  very  strongly  influenced  by  the  gas-drop  relative  velocity, 
detailed  information  on  the  aerodynamics  of  the  flow  in  a given 
combustor  configuration  is  essential  for  its  prediction.  Such 
information  can  be  obtained  either  from  measurements  or  from  simulation  of 
the  events  taking  place  in  actual  systems.  Owing  to  the  advances  made  in  the 
field  of  mathematically  modelling  turbulent  flow  phenomena  and  the  develop- 
ment of  powerful  techniques  for  solving  the  partial  differential  equations 
of  fluid  mechanics^^,  the  latter  approach  is  on  the  way  to  becoming  more 
efficient  and  economical. 


The  flow  equations  can  be  solved  using  cither  two  or  three  dimensional 

models^^,  however  the  solution  of  the  three  dimensional,  two  phase 

equations  require  considerable  computer  time  as  is  shown  in  Ref  12. 

In  practical  combustors,  the  region  in  the  general  vicinity  of  the 

fuel  injector  is  of  particular  concern, since  this  is  the  region  where 

soot  may  be  formed  from  unbumed  fuel  if  the  mixtures  are  locally  rich. 

(The  mechanism  of  soot  formation  is  discussed  in  Ref.  13  where  equations 

are  given  for  both  the  first  stage  representing  the  formation  of 

radical  nuclei,  and  the  second  stage  representing  soot  particle  formation 

from  these  nuclei).  Fig.  2 shows  the  flow  pattern  in  the  primary 

zone  of  a simple  gas  turbine  combustor,  computed  by  3-dimensional 
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mathematical  modelling  and  it  can  be  seen  that  the  flow  is  effectively 
2-dimensional  in  the  region  of  the  fuel  injector.  It  will  be  shown 
below  that  the  fuel  trajectory  lies  in  this  region  and  hence  2-dimensional 
modelling  of  the  spray  evaporation  is.  justified. 

The  main  purpose  of  the  present  study  is  to  demonstrate  a technique 
for  predicting  the  rate  of  evaporation  of  the  fuel  spray  by  combining 
the  abovementioned  methods  of  handling  turbulent  flow  situations  with 
spray  dynamics.  The  basic  approach  is  to  determine  the  flowfield  first, 
and  then  proceed  to  the  calculation  of  the  trajectories  of  individual 
droplets,  the  droplet  velocities  along  these  trajectories,  the  change  of 
size  distribution  and  finally  the  rate  of  evaporation. 

Separate  calculations  (to  be  published)  have  shown  that  at  conditions 
in  typical  gas  turbine  and  industrial  combustors,  the  influence  of  the 
droplets  on  the  flowfield  is  negligable  except  in  the  initial  region  of 
spray  break-up.  This  interaction  only  becomes  significant  when  the 
droplet  velocities  are  an  order  of  magnitude  greater  than  the  gas  velocity. 


ft  Sr 

is  I 


This  condition  sometimes  arises  when  fuel  is  injected  into  a recircu- 
lating flow  of  particularly  low  velocity,  but  is  unimportant  for  the 
cases  discussed  herin.  The  current  trend  towards  air  blast  atomizers 
and  other  devices  producing  fine  sprays  itrp lies  that  the  interaction  will 
continue  to  be  negligable  for  most  practical  combustors. 

To  demonstrate  the  method,  the  calculation  procedure  is  applied 
to  the  case  of  injection  into  a recirculation  zone  behind  a circular 
baffle  in  a duct.  The  governing  equations  of  the  fluid  flow  together 
with  the  turbulence  model  employed  and  the  method  of  solution  are 
described  in  Section  2.  Section  3 deals  with  some  aspects  of  spray 
evaluation  and  the  results  are  presented  in  Section  4.  In  Section  5 
the  results  of  various  experimental  tests  are  interpreted  in  terms  of 
the  theoretical  predictions. 


4 - 





2.  The  Theoretical  Formulation 

In  the  absence  of  heat  and  mass  transfer  and  chemical  reactions 
the  governing  equations  of  the  flow  are  the  well  established  equations 
for  the  conservation  of  mass  and  momentum,  supplemented  by  the  equations 
of  the  kinetic  energy  of  turbulence  and  the  rate  of  dissipation  of 
energy,  which  constitute  a simultaneous  nonlinear  set.  In  the  case  of 
two  dimensional  and  axially  symmetric  flow  problems  it  is  advantageous 
to  work  with  vorticity  and  stream  function  as  the  dependent  variables 
rather  than  the  components  of  velocity  and  pressure.  The  advantages 
exhibit  themselves  in  the  elimination  of  pressure  from  the  equations 
and  the  reduction  of  the  total  number  of  equations  to  be  solved  from 
five  to  four.  The  primitive  variables  are  recovered  after  the  solution 
of  the  problem  is  established. 

The  general  form  of  the  equations  in  cylindrical  coordinates, 
as  given  by  Gosman  et.  al.^^  is  the  following: 


a.  - 

41  ox  9y 


Here,  ^ is  the  general  dependent  variable  which  can  stand  for 
vorticity  divided  by  the  distance  from  the  axis  of  symmetry  lo/y,  the 
stream  function  \l),  the  kinetic  energy  of  turbulence  k,  and  the  rate  of 
energy  dissipation  e.  The  corresponding  expressions  for  the  quantities 
a^,  b^,  c^  and  d^  are  given  in  Table  1. 
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Table  1.  Expressions  for  quantities  appearing  in  eq.  (1). 
The  effective  viscosity  is  the  sum  of  laminar  viscosity  and  a 
turbulent  counterpart  originating  from  the  mixing  length  hypothesis. 
The  turbulent  viscosity  is  linked  to  the  kinetic  energy  of  turbulence 
and  the  energy  dissipation  rate  in  the  following  fashion: 


(2) 


Here,  p is  the  density  of  the  fluid  and  C , C„,  C- , C-  a,  and 

p D 1 z K,err 

0^  are  constants,  the  recommended  values  for  which  are  given 
in  Table  2. 
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Table  2.  Recommended  values  of  the  constants  of  the 
turbulence  model  (14) . 
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The  details  of  placing  the  general  partial  differential  equation 
(1)  into  finite  difference  form  by  integrating  it  over* a control  volume 
is  given  by  Gosman  et.al.^^.  A typical  control  volume  and  associated 
symbols  are  shovm  in  Figure  3.  The  successive  substitution  formula 
which  can  stand  for  either  of  the  dependent  variables  is  of  the 
following  form: 


K ' Ve  * Vw  * 


where  the  C*s  are  the  coefficients  which  contain  contributions  from 
the  convective  and  diffusive  fluxes  and  S is  the  source  term. 

Due  to  the  elliptic  nature  of  the  equations,  the  successive  sub- 
stitution formula  must  be  supplemented  by  boundary  conditions 
specified  at  all  boundaries.  Attention  shall  now  be  directed 
towards  the  specification  of  boundary  conditions  relevant  to  the 
baffle  problem  at  hand. 

The  geometry  of  the  domain  of  integration  is  shown  in  Figure 
4.  Four  different  types  of  boundary  can  be  distinguished:  the  inlet, 
the  walls,  the  axis  of  symmetry  and  the  exit.  The  conditions  at  these 
boundaries  will  be  handled  separately. 

InJizJt.  The  conditions  at  the  inlet  are  specified  once  and  for  all. 

In  other  words  the  values  of  the  dependent  variables  are  fixed  there 
and  do  not  require  updating  during  the  Iteration  process. 

^oJUi  Bou.nda/u.t&.  The  stream  function  does  not  present  difficulties 
and  its  value  is  constant  along  the  walls.  However,  it  is  known  from 
practice  that  steep  gradients  of  transport  properties  and  other 
dependent  variables  occur  near  the  walls  in  turbulent  flow.  To 
avoid  fine  grids  and  excessive  computing  time  use  is  made  of  the 


wall  flux  relations  near  the  walls.  The  familiar  log-law  of  the 
wall  can  be  stated  as: 


(4) 


where,  vt  is  the  laminar  viscosity,  p is  the  fluid  density,  n is  the 

normal  distance  from  the  wall,  u is  the  streamwise  velocity  at  a 

P 

point  near  the  wall,  u^  is  the  'shear  velocity*  and,  E and  k are 
constants  whose  values  are  9.793  and  0.4187  respectively.  The 
boundary  conditions  for  all. the  dependent  variables  except  the 
stream  function  are  derived  from  equation  C4  ) with  the  assumption 
that  the  flow  is  essentially  one  dimensional  in  the  vicinity  of 
the  walls.  These  conditions  are  specified  at  next  to  wall  grid 
points. 

AxaM  ol  SytmeJJuj.  The  stream  function  must  be  constant  along  the 
axis  of  symmetry  and  its  value  is  usually  taken  as  zero.  Also,  it 
is  obvious  that  the  normal  derivatives  of  k and  t must  vanish  on  the 
axis.  On  the  other  hand,  although  the  vorticity  itself  is  zero  the 
same  is  not  valid  for  u/y.  A satisfactory  boundary  condition  for 
this  variable  can  be  derived  by  expanding  into  a power  series  in 
y for  constant  x in  the  immediate  neighbourhood  of  the  axis  and 
substituting  this  into  the  definition  of  vorticity.  This  treatment 
leads  to  the  conclusion  that  the  normal  derivative  of  m/y  is  also 
zero  at  y equals  zero. 

ExaX.  Since  the  distributions  of  the  dependent  variables  at  the 
exit  plane  are  not  known  beforehand,  it  is  necessary  to  employ 
gradient  type  boundary  conditions  there.  This  is  expressed  mathe- 
matically as  3^/3x  = 0. 
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The  set  of  simultaneous  linear  equations  together  with  the 

algebraic  relationships  which  define  the  conditions  at  the  boundaries 

are  solved  using  the  Gauss-Seidel  iterative  scheme.  To  suppress 

divergence  provoked  by  large  variations  in  the  source  terms  of  the 

successive  substitution  formula  (.3  ) , k,  e and  p are  underrelaxed. 

The  underrelaxation  factors  for  tne  kinetic  energy  of  turbulence  k 

and  the  energy  dissipicion  rate  e are  set  at  0.7,  and  that  of  turbulent 

viscosity  at  0.5.  Other  variables  do  not  require  underrelaxation-. 

The  iteration  process  is  terminated  when  the  maximum  percentage 

change  of  all  dependent  variables  in  successive  iterations  falls 
-4 

below  10  . 


3.  Spray  Evaluation 

Once  the  flow  field  is  established  by  the  methods  described 
in  the  previous  sections,  the  next  step  is  the  determine  the  time 
rate  of  evaporation  of  a fuel  spray  into  this  field  at-  a given 
point.  (N.B.  Alternatively  the  flowfield  could  be  determined 
experimentally) . 

The  equation  of  motion  of  a droplet  in  the  absence  of  all 
external  forces  except  drag  is: 

dt  D p ' P ' (5) 

where  and  u are  the  absolute  droplet  and  gas  velocities  res- 
pectively, t is  the  time,  is  the  liquid  density,  is  the  gas 
density,  is  the  drag  coefficient  and  D is  the  droplet  diameter. 
When  written  in  component  form  this  becomes, 

du  _ 

^ P^B  " (» 

for  the  x-direction,  and 

dv  _ 

dt'’  = ■ ¥ -^2  < V 

f'l”  (7) 


for  the  y-direction.  Here  p is  the  viscosity  of  the  gas  and  Re  is 

the  relative  Reynolds  number  defined  as  follows: 

Re  = Dp  lu  - uj /p  , 
g'  P ' 

If  the  droplet  is  evaporating  a third  equation  giving  the  rate  of 
change  of  its  diameter  with  respect  to  time  is  needed.  In  the  case 
of  forced  convective  evaporation  this  can  be  written  as, 

^ ( 1 + 0.23  Re"'^^)  . /gN 

dt  2D 


is  the  evaporation  rate  constant  whose  value  depends  on  the  type 
of  fuel  as  well  as  the  properties  and  temperature  of  the  surrounding 
medium.  A widely  used  expression  for  derived  from  the  quasi-steady 
analysis  of  droplet  combustion  is  that  given  by  Wise  and  Agoston^^: 

p 8A  . , 

C = log  {l  + — c(t  - T H 

b P^Cp  V (9) 


where,  A and  Cp  a:re  the  thermal  conductivity  and  the  specific  heat 
at  constant  pressure  of  the  surrounding  gas,  L is  the  latent  heat 
of  evaporation,  and  are  the  ambient  and  drop  surface  teii?>eratures 
respectively.  This  equation  shows  that  the  evaporation  constant  is 
strongly  influenced  by  the  properties  of  the  gas,  A/Cp,  and  that  it 
depends  logarithmically  on  the  heat  of  vapourization  and  temperature. 

Equations  (6),  (7)  and  (8)  constitute  a set  of  simultaneous 
ordinary  differentials  equations  which  are  supplemented  with  the 
auxilliary  relations  for  the  evaporation  rate  constant  and  the  drag 
coefficient,  and  subject  to  a given  set  of  initial  conditions.  The 
components  of  gas  velocity  needed  for  the  solution  of  these  equations 
are  recovered  from  the  distribution  of  the  stream  function  and  are 
specified  at  the  grid  nodes.  Since  the  distrubutions  of  these  are 
not  available  as  continuous  function  of  x and  y,  it  is  assumed  that 
u and  V are  uniform  within  the  cells  and  that  their  values  are  the 
same  as  those  prevailing  at  point  P which  the  cell  encloses.  An 
alternative  but  probably  time  consuming  approach  would  be  to  express 
u and  V as  analytical  functions  by  a least  squares  fit  or  any  other 
similar  process. 

Focusing  attention  on  the  translation  of  a droplet  in  a given 

cell  and  considering  only  a very  short  segment  of  its  trajectory, 

1/2 

it  can  be  stipulated  that  the  quantities  Cj^Re  and  (1  + 0.23  Re  ) 


vary  only  negligibly  throughout  that  segment.  It  is  then  possible 
to  integrate  equations  (6),  (7)  and  (8)  analytically.  However,  it 
must  be  remembered  that  the  solutions  thus  obtained  are  only  valid 
i£  the  time  increment  under  consideration  is  sufficiently  small. 
Integration  of  equation  (8)  gives, 

- C,  (1  + 0.23  Re^^^)(t  - t ) 

0 D u 


for  the  time  dependence  of  the  droplet  diameter.  Re  is  the  relative 
Reynolds  number  evaluated  at  time  t^.  Substituting  this  into  equations 
(6)  and  (7)  and  integrating. 


V = V + (v  - v)(D/D.)' 
P P»0  0 


“ " (12) 

are  obtained  for  the  components  of  the  droplet  velocity.  The  co- 
ordinates of  the  instantaneous  location  of  the  drop  are  recovered  by 
integrating  these  with  respect  to  time.  Giving, 

Ti  2((i+1 ) 


I’n  r,  2(a+l) 


6 * C (I  + 0.23  . a)  . 


The  set  of  equations  (11) > (12),  (13)  and  (14),  which  give 
the  velocity  and  the  location  of  a droplet  of  diameter  D at  a given 
instant  of  time,  must  be  supplied  with  the  relationship  between  the 
drag  coefficient  and  the  relative  Reynolds  number.  The  three  point 
formula  of  Morsi  and  Alexander^^  has  been  eiiq>loyed  for  the  evaluation 
of  the  drag  coefficient.  This  formula  is  expressed  as. 


^1  ^2 

C = — + ^ + K 
Re  Re" 

where  K^,  and  are  constants  the  values  of  which  depend  on 
the  working  range  of  Reynolds  number. 

Since  the  spray  is  made  up  of  a spectrum  of  droplet  sizes, 
information  on  the  initial  size  distribution  is  also  required  for 
the  estimation  of  the  rate  of  its  evaporation.  This  information  is 
normally  obtained  from  measurements.  An  extremely  rapid  and  simple 
method  which  can  be  used  for  determining  the  size  distribution  in 
terms  of  the  Rosin-Rammler  parameters  is  the  laser  diffraction 
technique  of  Swithenbank  et.  al.^^.  Owing  to  its  simplicity  the 
Rosin-Rammler  equation  has  been  employed  in  this  study,  however 
other  functions  such  as  the  Nukiyama-Tanasawa  or  the  log-normal 
distributions  can  be  accomodated  when  required. 

The  calculation  procedure  is  further  facilitated  if  the  dis- 
tribution curve  is  divided  into  an  adequate  number  of  size  groups 
each  represented  by  a mean  diameter  D.  The  initial  mass  fraction 
of  each  group  can  then  be  expressed  as. 


Am  = exp{-(D/^/}  AD 

^ (16) 

where  Am  is  the  mass  of  drops  of  size  ±LD/2  around  D per  unit 

volume  of  injected  fuel,  D is  the  size  constant  and  5 is  the 

distribution  parameter.  Normally  15  to  20  size  intervals  are 

required  to  represent  the  spray  adequately.  For  each  size  range 

the  following  steps  are  performed  at  every  time  increment: 

1.  The  relative  Reynolds  number  is  calculated  using  the 
conditions  prevailing  at  the  beginning  of  the  increment  denoted 
by  time  t. 

2.  The  drag  coefficient  is  obtained  from  the  three  point 

1/2 

formula  and  the  quantities  Cjjffe/24  and  (1  + 0.23  Re  ) are 
formed. 

3.  The  drop  diameter  at  time  t+At  is  calculated  from  equation  10. 

4.  The  mass  remaining  at  t+At  is  found  from 

5.  The  new  values  of  u and  v are  calculated  from  equations  li. 
These  are  required  to  determine  the  quantities  involving  the 
Reynolds  number  for  the  next  increment  of  time. 

6.  The  location  of  the  droplet  at  time  t+At  is  obtained  from 


equations  14.  This  location  is  checked  against  the  cell 
boundaries  to  find  out  whether  the  drop  has  moved  into  one 
of  the  neighbouring  cells.  If  so  the  corresponding  values 
of  the  components  of  air  velocity  are  substituted  in  equations 
13  and  14. 


7.  Steps  1 to  6 sre  repeated  entil  the  drop  diameter  falls  to  zero 
or  impingement  on  wall  boundaries  occurs. 

Once  the  evaporation  histories  of  individual  size  groups  are  cal- 
culated the  mass  fractions  are  sinmmd  at  predetermined  instances  to 
obtain  the  total  evaporation  tine  curves.  It  can  be  noted  that  the 
above  formulation  also  yields  drop  trajectories,  the  drop  velocity 

histories  along  these  trajectories  and  the  size  distribution  curves 
different  times. 


4.  Results 

The  calculations  are  performed  for  the  geometry  of  an  existing 
model  combustor  can  of  internal  diameter  8 cms.  and  length  24  cms.. 

The  baffle  which  is  situated  at  the  inlet  is  hemispherical  in  shape, 
its  diameter  being  half  that  of  the  can.  The  air  inlet  velocity  was 
taken  as  10  m/sec. 

Figurue  5.  shows  the  streamlines  of  the  mean  flow  in  the  recir- 
culation zone  and  the  development  of  the  axial  velocity  distribution. 

It  can  be  seen  that  the  bubble  extends  to  about  one  baffle  diameter. 

The  inspection  of  the  velocity  profiles  reveals  that  the  velocities 
in  the  recirculation  zone  are  an  order  of  magnitude  smaller  than  those 
in  the  annular  jet.  Generally,  these  results  agree  at  least  qualitatively 
with  experiment.  The  distribution  of  the  kinetic  energy  of  turbulence 
in  the  wake  is  given  in  Figure  6.  The  conclusion  can  be  drawn  that  the 
level  of  turbulence  exhibits  a maximum  in  the  latter  part  of  the  bi'.bble. 

Again,  these  results  show  a qualitative  agreement  with  the  experimental 

18 

findings  of  Davies  and  Beer,  although  the  two  cases  are  not  entirely 
similar^ 

The  injected  spray  has  been  taken  to  have  a total  cone  angle  of 

45°.  All  the  droplets  were  assumed  to  leave  the  nozzle  with  an 

initial  velocity  of  21  m/sec.  The  calculations  were  performed  for 

two  different  size  constants  of  the  Rosin-Rammler  distribution,  50 

and  100  microns,  and  the  distribution  parameter  was  set  at  2. 

The  values  of  the  evaporation  constant  worked  out  from  equation  (9) 
o “A 

for  T = 1500  K and  A/c  = 7.8  x 10  grms./cm-sec  range  from  0.002 

00  p 

2 

to  0.01  cm  /sec  depending  on  the  properties  of  the  fuel.  Since  the 


effect  of  the  oxidizer  mass  fraction  is  omitted  in  equation  (9) , 


these  are  lower  than  the  corresponding  burning  rate  constants.  An 

2 

overall  value  of  0.005  cm  /sec  has  been  used  throughout  the  present 
calculations. 

The  trajectories  of  individual  droplets  are  presented  in  Figure 
7.  These  show  a slight  deviation  from  the  nominal  cone  in  the 
recirculation  region  and  bend  thereafter  towards  the  downstream 
direction.  It  can  be  seen  that  the  smaller  droplets  evaporate 
coii^)letely  in  the  vicinity  of  the  point  of  issue,  whereas  those  with 
initial  diameters  larger  than  110  microns  hit  the  wall  of  the 
duct.  Figures  8 and  9 show  the  profiles  of  unevaporated  mass 
fraction  against  time  for  two  different  values  of  the  size  constant. 

As  expected,  forced  convection  has  an  appreciable  influence  on  the 
rate  of  evaporation.  Again,  both  figures  show  that  this  rate  increases 
with  decreasing  size  constant.  Figure  10  gives  the  change  of  the 
distribution  function  with  time  for  D = 100  microns.  The  peak  of 
the  curve  shifts  to  the  right  as  time  elapses. 


4.1  Heat  Transfer  to  Droplets 

The  evaporation  lag  of  droplets  which  is  defined  as  the  time  required 
to  bring  the  surface  temperature  of  a droplet  to  the  boiling  point  of 
liquid  is  considered.  The  time  rate  of  change  of  the  droplets  temperature 
is  given  by 


6 Nu  k 

n2 

p,D  c , 

1 P,d 


(T  " 
00 


uhere>  k is  the  thermal  conductivity  of  the  gas,  Nu  the  Nusselt  number, 

Cp  j the  specific  heat  of  droplets  and  T^  the  temperature  of  the  surrounding 
medium.  Under  forced  convective  motion  the  Nusselt  number  is  given  by  the 
following  equation  which  is  due  to  Ranz  and  Marshall, 

Nu  = 2 + 0.6  Re  Pr 

where,  BE  is  the  relative  Reynolds  number  and  ?r  is  the  Prandtl  number. 
Supposing  T is  constant  throughout,  the  above  equation  can  be  integrated 


to  give 


T = T^  - (T„  - T^)  ^ *^'^'^1°  '^p,d^ 


where,  is  the  initial  droplet  temperature  and  D is  its  diameter. 

A calculation  has  been  performed  with  T = 1500°K,  T^  = 320°K, 

k * 10  ^ W/cm  °K,  c j = 3 Ws/gr  °K  the  results  of  which  are  presented 
p,a 

in  Figs  11  & 12.  The  boiling  point  of  fuel  was  taken  to  be  450 °K. 


1. 


5.  Discussion. 

The  results  presented  above  predict  several  significant  features. 

In  the  hot  combustion  environment  sprays  of  the  type  commonly  used  in 
combustors  tend  to  travel  in  near  straight  lines.  The  smaller  droplets 
< 50  pm  typically  only  travel  2-3  centimetres  before  evaporating 
away,  whereas  droplets  greater  than  100  ym  travel  many  centimetres  and 
may  hit  the  combustor  wall.  In  general,  as  a droplet  becomes  small 
enough  to  be  deflected  by  the  gas  flow,  its  life  is  so  short  that  it 
quickly  vanishes.  This  conclusion  is  also  true  when  allowance  is  made 
for  the  heating-up  time  of  the  drop  as  this  factor  only  contributes 
about  10%  to  the  range.  These  results  account  for  the  experimental 
observations  at  Sheffield  University,  (Chigier  et  al  ) in  which 
very  detailed  measurements  were  made  of  the  streamlines,  droplet 
'trajectories',  temperatures  and  concentrations  in  a flame  stabilised 
behind  a disc  furnished  with  a hollow  cone  spray  nozzle  of  41°  total  cone 
angle.  The  overall  spray  model  deduced  in  Ref  19  is  reproduced  in 
Fig.  13.  In  the  photographic  method  used,  it  is  difficult  to  follow 
the  time  history  of  individual  droplets  since  their  size  changes 
continuously  as  they  evaporate.  Nevertheless  the  experimental  results 
clearly  show  that  in  the  hot  environment,  droplets  larger  than  100  ym 
move  in  straight  lines,  whilst  as  they  near  the  end  of  their  life  they  tend 
to  align  with  the  local  streamlines.  In  the  cold  unbuming  environment 
the  drag  on  the  particles  is  much  larger  and  their  lifetimes  are  much 
longer.  Consequently,  small  droplets  may  be  captured  in  the  turbulent 
recirculation  zone  and  are  then  observed  moving  in  all  directions  in 
this  region. 
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Flame  stabilization  behind  a diae  in  a duct  using  a 90°  cone 
angle  liquid  propane  fuel  spray  has  been  studied  at  Perdud  Oniveraity“. 
Detailed  concentration  and  temperature  maps  are  presented  from  which 
the  spray  path  and  reaction  zones  are  deduced.  Again  the  results 
agree  with  the  predictions  of  the  analysis  presented  above,  and  the 
relationships  estahlished  show  how  the  evaporation  rate  can  be  related 
to  pollutant  formation  in  the  various  combustion  zones. 

Finally  the  fact  that  large  droplets  arrive  at  the  predicted  location 
on  the  walls  of  a gas  turbine  combustor  can  be  experimentally  demonstrated 
by  simply  inspecting  the  wall  of  the  burning  combustor  using  a very  high 

intensity  lamp  to  overcome  the  flame  radiation.  The  intermittently  wet 
surface  spots  may  be  readily  seen. 
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Measured  stability  loops  for  70  mm  diam.  gas  turbine 
combustor. 

Axial  and  radial  velocity  distributions  in  primary  zone  of 
70  mm  diameter  combustor. 

a)  In  plane  10  mm  downstream  of  fuel  nozzle 

b)  In  plane  6 mm  upstream  of  primary  jets. 

A Typical  volume  element. 

The  domain  of  integration  and  boundary  conditions 

The  streamlines  of  the  mean  flow  in  the  recirculation  zone 
and  axial  velocity  distributions  at  several  downstream 
locations. 

Lines  of  equal  turbulence  kinetic  energy  in  the  wake. 

Trajectories  of  individual  droplets.  Arrows  indicate  burnout 
locations.  Solid  line  is  the  nominal  cone. 

Plots  of  unevaporated  mass  fraction  against  time,  without 
forced  convection  (fi  =2). 

Plots  of  unevaporated  mass  fraction  against  time,  including 
forced  convection  (6  « 2). 

The  change  of  distribution  function  for  size  constant,  D = 100 
microns  and  distribution  parameter  6 = 2 at  t = 0. 

Rate  of  evaporation  with  and  without  allowance  for  droplet 
heating  up  time. 


Droplet  trajectories  with  and  without  allowance  for  droplet 
heating-up  time. 

Physical  model  of  spray  burning  in  the  wake  of  a stabilizer 
disc^^. 


Axis  of 
symmetry 


Fig.  3.  A typical  volume  element. 


~ '^top  * values  of  w/y,  k and  e 

from  wall  fiuictions 


8 4.  The  domain  of  integration  and  boundary  conditions 


axial  velocity  distributions  at  several  downstream 


Plots  of  unevaporated  mass  fraction  against  time,  without 
forced  convection  (6  «*  2) . 


Plots  of  unevaporated  mass  fraction  against  time,  including 
forced  convection  (6  * 2) 
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